Downfolded Self-Energy of Many-Electron Systems 
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Starting from the full many-body Hamiltonian of interacting electrons the effective self-energy 
acting on electrons residing in a subspace of the full Hilbert space is derived. This subspace may 
correspond to, for example, partially filled narrow bands, which often characterize strongly cor- 
related materials. The formalism delivers naturally the frequency-dependent effective interaction 
(the Hubbard U) and provides a general framework for constructing theoretical models based on 
the Green function language. It also furnishes a general scheme for first-principles calculations of 
complex systems in which the main correlation effects are concentrated on a small subspace of the 
full Hilbert space. 
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One of the most important quantities in studying the 
electronic structure of materials is the spectral function, a 
quantity that contains information about the electronic 
excitations of the system obtained from photoemission 
and inverse photoemission experiment corresponding to 
the addition and removal of an electron. A suitable 
tool to describe the spectral function is the one-electron 
Green's function. To calculate the Green function of the 
full many-body Hamiltonian is a tremendous task and 
only recently such calculations, within the GW approx- 
imation (GWA) [H, 0, are feasible thanks to the rapid 
progress in computer performance. However, two seri- 
ous difficulties are hampering further progress: firstly, for 
many cases, the system size is too large for realistic GW 
calculations and secondly, from a theoretical point of 
view, the GWA may not be sufficient to treat correlation 
problems, especially those in materials with strong cor- 
relations containing partially filled narrow bands. These 
materials with many intriguing properties have been dis- 
covered in recent years Q. 

The traditional approach for treating such systems is 
to introduce a model Hamiltonian, focusing on a small 
subspace of the full Hilbert space that is considered to 
be most relevant for the correlation problem at hand. No- 
table examples are the Hubbard model and the Anderson 
impurity model with the assumption of a local statically 
screened Coulomb interaction known as the Hubbard U, 
which is treated as an adjustable parameter of the mod- 
els. These models have given us a lot of physical insight 
into correlation problems ranging from the Kondo effect 
to high temperature superconductivity. The presence of 
an adjustable parameter U, however, limits the predic- 
tive power of the models, preventing accurate and pre- 
cise quantitative calculations. Moreover, these models 
are not strictly derived from the many-body Hamilto- 
nian, but rather they are intuitively postulated, which is 
not satisfying from theoretical point of view. 



The purpose of the present work is to derive from the 
full many-body Hamiltonian an exact expression for the 
self-energy corresponding to a subspace of the full Hilbert 
space, expressed in terms of the one-particle Green's 
function of the subspace. The motivation for this is that 
in many physical problems it is too complicated to treat 
the full Hilbert space on equal footing. A sensible and 
physically motivated approach is to identify a subspace, 
where most of the correlations take place. This subspace 
could be identified as, for example, partially filled nar- 
row bands crossing the Fermi level, a typical situation in 
materials with strong correlations. The task is then to 
derive an effective self-energy acting on electrons resid- 
ing in this subspace. As a result of the derivation, the 
effective screened interaction among electrons belonging 
to the subspace, i.e., the Hubbard U, emerges naturally. 

Let us start by defining some basic variables and quan- 
tities. The many-electron Hamiltonian is given by 



H = J d 3 r ip + (r)h (r)ip(r) 

+ i J d 3 rd 3 r' V + (r)V> + (r>(r - r')^(r')^(r), (l) 

where ho is the one-particle part of the Hamiltonian and 
is the field operator. We aim at downfolding the 
many-body problem onto a subspace of the full Hilbert 
space, which can consist of, for example, 3d or 4f orbitals. 
From now on we refer to the subspace as the d subspace. 
We first divide the complete field operator into the d field 
and the rest, denoted by d and r, respectively: 



V>(r) = ip d (r) + Vv(r) 



d 



Xd{r)c d + ^2xAr)c r . (2) 



Xd and Xr are the one-particle orbitals and Cd and cy 
are the associated annihilation operators. We use a con- 
vention that (rt) is represented by a number. Since we 
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will utilize the Schwinger functional derivative technique 
[i| to develop a closed set of equations for G d , where a 
probing field ip(l) is applied to probe the linear response 
of the Green function, we define the following Green's 
functions in the Dirac or interaction representation 



are summed and repeated variables are integrated, unless 
they appear on both sides of the equation. From the 
Heisenberg equation of motion for the field operators ijjd 
and we obtain the equations of motion for G d and 

nrd. 



iG d (l,2) = 



iG rd (l,2) 



(T [5^(1)^(2)]) 
(S) 

(S) 



S = T exp 



dt / d A r tp(rt)ip + (rt)?p(rt) 



(3) 



(4) 



(5) 



Before proceeding further with a detailed derivation, 
we first summarize schematically our final result of a 
closed set of equations for the effective self-energy and 
the Green function of the d subspace. 



£ = S d + £ rd + Z drd , 
P d = -iG d TG d , 

r=1 + ^ GdrGd ' 
W = W r + W r P d W 

G d = g d + g d Y>G d . 



(6) 
(7) 

(8) 

(9) 
(10) 



In place of the bare Coulomb interaction v, we have in © 
W r , a frequency-dependent effective interaction among 
electrons living in the chosen d subspace. This effective 
interaction is screened by P , the polarization propaga- 
tor of the d electrons with a vertex T, yielding the fully 
screened interaction W of the full system. T is an ef- 
fective vertex acting on the d subspace and consequently 
different from the full vertex. If the subspace is isolated 
and identified as the Hilbert space of a Hubbard model 
and W r is approximated by a local and static interac- 
tion (the Hubbard U), S d is then the exact self-energy 
of the Hubbard model. The explicit expressions for the 
self-energies will be shown later. In addition to S d the 
true effective self-energy £ contains a term, E rd , which 
is the self-energy arising from the r subspace acting on 
the d subspace and a hopping term £ drd representing hy- 
bridization between the d and r subspaces and scattering 
processes of an electron from the d to the r and back to 
the d subspace. Thus, E is not merely a projection of 
the full self-energy onto the d subspace. As will be ex- 
plained later, for a given W r the effective self-energy is 
a functional of G d only. The last equation is the Dyson 
equation with a noninteracting Green's function g d of the 
d subspace. 

We now proceed with the derivation of the above set of 
equations and use the convention that repeated indices 



i J-G d (l, 2) - A d (l, 3)h°(3)[G d (3, 2) + G rd (3, 2)] 

Oti 

- A d (l, 3) [M d (3, 4) + M rd (3, 4)] G d (4, 2) = A d (l, 2), 

(11) 



i — G rd (l, 2) - A r (l, 3)h°(3)[G d (3, 2) + G rd (3, 2)] 

Oti 

- A r (l, 3) [M d (3, 4) + M rd (3, 4)] G d (4, 2) = 0, (12) 
where 

A i (l,2) = Xi(ri)x*(r 2 )5(t 1 -t 2 ), (i = d,r), (13) 
and we have defined the mass operators according to 

A/ d (l,4)G d (4,2) 

= -iv(l - 3) (T [^+(3)^(3)^(1)^(2)] ) , (14) 



M rd (l,4)G d (4,2) 

= -iv(l - 3) (T [V- + (3)V(3)^(l)V + (2)]) . 



(15) 



It is important to realize that the mass operator M rd is 
defined with respect to G d , rather than G rd . Noting that 
SS/5ip(3) = -iT [5ty>+(3)V>(3)] the four-field term can be 
related to SG d /S(p and 6G rd /5ip [ESS]. Thus - de&n ~ 
ing the self-energy as the mass operator less the Hartree 
potential Vjj we obtain 



S^(l,4)G d (4,2) = ? V^(3,l) ^ ( ( 1 3) 2) , (j=d,rd) (16) 

We have defined V = Vh + <p and the screened Coulomb 
interaction W — (5V/Sip)v — e~ 1 v. 

To derive an effective self-energy operator for G d , we 
need to eliminate G rd . At this stage it is convenient to 
work in the basis representation. We denote a basis in 
the d subspace by Greek letters (a) and in the rest of 
the subspace by small roman letters (n) whereas capital 
letters (N) denote the full Hilbert space. Note also that 
G d has only matrix elements in the d subspace whereas 
G rd has only off-diagonal elements between the d-r sub- 
spaces. However, the self-energies E d and T, rd can have 
matrix elements between the d-d and d-r subspaces. 

The equations of motion of G d and G rd in frequency 
space become, 
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uS aj -h aj -^ d ai (Lo)-^(u)]G^(io)-h an G%(cu) = S aP , 

(17) 



,,„, -/ lnm ]G^( W )-K 7 +Sf l7 (c)+S^(c)]G^(c) = 0, 

(18) 

where the Hartree potential has been absorbed into h — 
h° + Vh- Note that h is static unless a probing time- 
dependent field <p(rt) is applied. Using the second equa- 
tion we eliminate G rd : 



= 9 r nm {u)[h mi + E^(w) + X^MIG^M, 

(19) 

where g r nm = [u> — h] n ^ n , which is a noninteracting Green's 
function of the r subspace. A small letter g is used to 
denote a noninteracting Green's function. The effective 
self-energy for G d is defined as follows. 



[u)5 al — h ai 



E^HIG^H = S af3 , (20) 



E Q » = E^» + E£( w ) + E^H, (21) 



Ej» = h an g r nm (u)[h mi + Sl» + E*(u;)]. (22) 

We are now in the position to calculate E d and Ti rd . 
We first note that G d {G d y 1 = A d since G d is only 
defined within the d subspace and 8Vm n {t) / 8V (r' t') = 
X* M (r')xN(r')S{t - t'). Multiplying both sides of OH) on 
the right by (G d )~ 1 , and using the identity G d S(G d )~ 1 + 
5G d (G d )- 1 = we obtain 



^%i3(ti~t 2 ) — iWN v .^v{h—tl)G d K (ti-t i )T K Jl(t4 : —t2',tz). 

(23) 

The vertex Y is given by 
Tf v (h-t 2 ,t 3 ) 

= S ail Sp v S{tl - t 2 )S(ti - t 3 ) 

H ^G d ^ti^ GdyK<yti _ <5 ^ r ^^ 5 _ i 6,^3)G^ I) (i 6 )- 

(24) 

The last line has been obtained by using a chain rule 
6E/SV = {5Y, d /8G d ){5G d /5V) and the identity SG d = 
-G d 5{G d )- 1 G d . 
Similarly 



E^^i — t 2 ) — iWKm,MN(h — h) 

x y-S ^-fa). (25) 

To calculate E rd we need to take the functional deriva- 
tive of G^ a with respect to the probing field. However, 
G^ a cannot be inverted and consequently the usual trick 
of using the identity S(GG~ 1 ) = does not apply. We 
therefore must calculate 8G rd /5V directly. To do this we 
use G rd (uj) in 033) j Fourier transform back into time do- 
main and perform the functional derivative. Calculating 
5G rd /5V and substituting the result in ([25]) yields 

= X%Y (r \h-t 3 )gl n (t 3 -t 5 ) 

x [h n/3 s(t 5 - 1 2 ) + E^(t s - 1 2 ) + s;^(i s - 1 2 )\ 

+ E^ (r) (ii - t 2 ) 

+ iWN mtftu (t 3 - ti)g r mn (tx - t 5 ) 

~ 6X? lp (t 5 -t 2 ) 5Z r n d (t 5 -t 2 ) ' 
SV^ih) SV^(t 3 ) 

+ iWN mtftu (t 3 - ti)g r mn {ti ~ t 5 ) 
x [h nX 8(t 5 - i„) + E^ A (t 5 - t 6 ) + E^(i 5 - i 6 )] 

xG^(t 6 -t 4 )r^(t4-t 2 ;t3), 

(26) 



where 



(27) 



Finally, the polarization propagator of the d subspace is 
given by 



5V{2) SV{2) 



(28) 



•fa/a,fti/(*i — *z) — -*G„ 7 (ti-t 3 )r^(i3-i4,t 2 )G^ ;3 (<4-i ] f ). 

(29) 

Writing the full polarization of the system as a sum of 
the d- and r-subspace polarization, P = P d + P r , it is 
straightforward to verify that with 



W r 



vP r W r = [1 -vP r \ 



(30) 



Eq. © is equivalent to W — v + vPW. This completes 
our derivation of the closed set of equations for G d and 
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the effective self-energy, which is schematically given in 
Eqs. ©, ©, ©, ©, and (10]), and in detail in Eqs. 
(121]), ((221), ©, and (2ni)- This set of equations may 
be regarded as a set of downfoldcd Hedin's equations for 
a given W r . This quantity may be calculated, for exam- 
ple, within the constrained random-phase approximation 
(cRPA) scheme Q. For a given screened interaction W, 
the self-energies and hence the effective self-energy are 
functionals of G d only, as can be seen in Eqs. ([23]) . ([26]) . 
and (|22|) . The dependence on G rd has been eliminated 
and there is no explicit dependence on G r . It is impor- 
tant to note that we have defined the effective self-energy 
such that 



T} aj G*g — £ F Q, 



7^7/3 



^anGnf: 



(31) 



where £ F and Q are respectively the self-energy and the 
Green function of the full Hilbert space. The first term 
on the right-hand-side is confined to the d subspace and 
the second represents a coupling between the d and r 
subspaces. Thus, the self-energy T, a p is not a simple 
projection of the full self-energy onto the d subspace but 
it contains the effects of the off-diagonal matrix elements 
of £ F and Q between the d and r subspaces. However, G d 
is equivalent to the projection of Q onto the d subspace 
as can be seen in (3]). 

It is instructive to compute the effective self-energy in 
our formalism within the GWA, where the vertex correc- 
tions 5Y, d /8V and SY, rd /5V are neglected. 



GW(d), s . v Gf(r), 



+ <4„M [/i„ A + ££» + £;! m] 

x Gf 1 (uj)W am ^f3{uj) 



(32) 



where E^ (r) and E^ d are given in J27]) and j22j), and 

= (*/2^(^ W)PW^')- (33) 

The first two terms correspond to the projection of the 
full GW self-energy onto the d subspace. As discussed 
in the previous paragraph, our effective self-energy con- 
tains the effects of the off-diagonal self-energy between 
the d and r subspaces: The third and fourth terms 
with the square brackets are the self-energy contribu- 
tion from the r subspace on the d subspace. The term 

9mn(u) Kx + £^ A M + G%{u) may be inter- 

preted as a self-energy correction to G d arising from the 
off-diagonal elements of h, E d , and E rd . This together 
with W amn i3(u) forms & GW self-energy correction. The 
last term Y, drd is the hybridization term representing hop- 
pings of electrons from the d to the r subspace and back 



to the d subspace. Using the G W self-energy and keeping 
W fixed it becomes feasible to solve the vertex equation 
in (|24|) since the subspace of interest is usually much 
smaller than the full Hilbert space. The resulting ver- 
tex T can then be inserted into the self-energy in ([2"3"|) to 
obtain an improved self-energy beyond the GWA. 

The set of equations for the downfolded self-energy 
provides a general framework for constructing theoret- 
ical models in the Green function language. It offers an 
alternative to model Hamiltonians, which cannot read- 
ily take into account frequency-dependent interactions. 
Rather than first mapping the full Hamiltonian to a 
model Hamiltonian with a static U and then solving the 
model using the Green function technique, the present 
formalism allows for a direct route to the self-energy with 
a full frequency-dependent effective interaction. When 
W r is approximated by a local and static value it be- 
comes evident that the set of equations may be inter- 
preted as equivalent to the Hubbard model when E rd and 
Yi drd are neglected. Consequently, these two self-energy 
terms have no counterparts in the Hubbard model. Two 
ingredients are therefore missing in the Hubbard model, 
namely, the frequency dependence of U, which can be 
important || , and the self-energy effects arising from the 
r subspace. In the Anderson impurity model, the or- 
bitals of the impurity may be viewed as forming the d 
subspace and the influence of the r subspace is included 
as a hybridization term, which may be identified as the 
first term in (f22|) . However, the model does not include 
the self-energy effects E rd arising from the r subspace 
and as in the Hubbard model, the effective interaction is 
local and static. 

Finally let us consider some possible general applica- 
tions. The present formalism can be used to improve 
the Hubbard or the Anderson impurity model by provid- 
ing better parameters and by including the missing self- 
energy terms and the effects of the frequency-dependent 
W r within, for example, the GWA. The formalism also 
provides a route for numerical simplification of the GWA, 
by focusing on a certain subspace and treating the rest in 
an approximate, e.g., mean-field scheme. This is highly 
desirable since applications of the GWA to complex sys- 
tems are computationally demanding. Also, the formal- 
ism may facilitate a tractable way of going beyond the 
GWA fully from first-principles by inclusion of vertex cor- 
rections on the given subspace only. 

In conclusion, we have derived from the full many-body 
Hamiltonian a closed set of equations for the effective self- 
energy acting on a subspace of the full Hilbert space. The 
effective frequency-dependent interaction or the Hubbard 
U appears naturally in this formalism. Since the equa- 
tions are exact, they provide a general framework for 
handling complex systems in which the main correla- 
tion effects are concentrated on a certain subspace, such 
is the case in many correlated materials characterized 
by partially filled narrow bands. It offers an alterna- 
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tive to model Hamiltonians by constructing theoretical 
models in the Green function language allowing for a di- 
rect access to the self-energy and a natural inclusion of 
frequency-dependent interactions, usually not accounted 
for in conventional model Hamiltonians. 

We thank Silke Biermann for many fruitful discussions 
and suggestions. 
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